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I .  Introduction 

In  recent  years  the  study  of  failure  of  aircraft  fuel  cells  caused 
by  encounter  with  ballistic  threats  has  become  increasingly  important. 
Projectiles  that  penetrate  fluid  filled  cells  cause  damage  many  times 
greater  than  that  incurred  by  impact  with  an  empty  cell.  This  inter- 
action is  called  Hydraulic  Ram.  Hydraulic  ram  may  be  studied  by  iden- 
tifying several  phases  which  are  named  after  the  dominant  projectile 
energy  dissipation  mechanism.  These  phases  and  their  effects  on  fuel 
cell  wall  response  occur  during  different  time  scales.  The  phases 
usually  identified  are: 

1.  a  penetration  phase  where  the  projectile  punches  a  hole  in  the 
entry  wall 

2.  a  shock  phase  during  which  projectile  protusion  into  the  fluid 
drives  a  hemispherical  shock  wave  into  the  tank  liquid 

3.  a  drag  phase  during  which  projectile  energy  is  dissipated  by 
form  drag  and  a  cylindrical  vapor  filled  cavity  is  formed 

k.     a  cavity  oscillation  phase  which  results  from  the  growth  and 
eventual  collapse  of  the  drag  phase  cavity 

This  report  presents  the  results  of  experimental  studies  designed 
to  substantiate  current  fuel  cell  internal  and  wall  pressure  prediction 
models  for  the  shock  and  drag  phases  of  hydraulic  ram,  and  to  measure  the 
entry  wall  response  to  these  loadings. 


II.  Shock  Phase  Pressures 

12  3 
Several  investigators  '  '  have  shown  that  projectile  protusion 

into  a  fuel  cell  produces  a  hemispherical  shock  wave  which  is  centered 

at  the  impact  point.   Projectile  momentum  and  path  after  entry  have  no 

observable  effect  on  the  shock  shape  as  it  propagates  into  the  fluid. 

Since  impacts  with  obliquity  also  exhibit  this  characteristic  it  was 

concluded  that  shock  phase  characteristics  may  be  modeled  as  the  result 

k 
of  a  point  energy  source  release.   The  theory  of  Yurkovich  is  based  upon 

this  assumption.   The  theory  assumes  a  rigid  entry  wall  and  shock  radius 

that  is  proportional  to  a  constant  power  of  time.   Spark  shadowgraph 

12  3 

experiments  '  '  have  shown  that  blunt  projectiles  (spheres,  shells  with 

hemispherical  noses,  and  warhead  fragments  for  example)  exhibit  a  shock 
radius  that  is  proportional  to  the  0.8  power  of  time  and  that  projectiles 
with  ogive  nose  shape  produces  a  shock  radius  that  is  proportional  to  the 
0.9  power  of  time. 

The  flow  field  was  calculated  by  coupling  the  one- dimensional  Rankine- 
Hugoniot  equations  with  the  Tait  equation  of  state  for  isothermal  compres- 
sible liquids  to  determine  shock  front  jump  conditions  as  a  function  of 
the  shock  Mach  number.  Using  the  normal  strong  shock  assumptions,  the  radial 
distribution  of  the  flow  parameters  behind  the  front  was  then  determined 
by  assuming  a  power  law  density  distribution  and  integrating  the  continuity 
equation  to  determine  the  radial  velocity  distribution.   The  momentum 
equation  was  then  integrated  to  determine  the  pressure  distribution. 
Given  a  shock  phase  energy  release  E  the  shock  Mach  number  as  a  function 
of  radius  was  determined  from  energy  considerations. 

The  theory  of  Yurkovich  gives  no  information  as  to  the  value  of  the 
energy  release  during  the  shock  phase.   Figure  II-l  shows  a  qualitative 


view  of  projectile  velocity  decay  during  hydraulic  ram.   It  is  known 
that  in  most  instances  only  a  portion  of  the  impact  energy  E.  transferred 
during  shock  formation  since  projectiles  travel  many  inches  into  the  tank 
before  coming  to  rest.   In  the  shock  formation  region  projectile  velocity 
drops  sharply  as  a  result  of  shock  formation.  Velocity  decay  after  shock 
formation  is  due  to  projectile  form  drag.   If  the  projectile  is  considered 
to  be  a  piston  in  the  early  stages  of  shock  formation,  the  rate  of  change 
of  fluid  volume  displaced  would  determine  the  momentum  exchange  between 
the  projectile  and  the  fluid,  hence  the  shock  strength.  It  may  be  assumed 
that  a  characteristic  time  for  shock  formation  would  be  that  required 
for  the  projectile  to  penetrate  the  fluid  an  axial  distance  equal  to  that 
at  which  the  projectile  maximum  cross- sectional  area  occurs.  The  shock 
formation  characteristic  time  is  no  longer  than  a  few  microseconds  for 
most  cases  of  interest.  It  may  also  be  deduced  that  the  projectile  Mach 
number  with  respect  to  the  fluid  has  an  important  role  in  determining 
the  strength  of  the  shock  and  therefore  the  amount  of  energy  deposited 
during  the  shock  phase.  An  analytical  expression  for  the  shock  phase 
energy  release  as  a  function  of  projectile,  wall,  and  fluid  properties 
is  not  available  to  date.  Information  about  the  shock  formation  process 
could  be  calculated  numerically  using  computer  codes  similar  to  those  used 
to  study  projectile  penetration  dynamics  of  multiple  element  targets.  Use 
of  this  method,  however  would  require  a  systematic  parameter  study  which 
would  involve  large  expenditures  of  computer  time.  Experimental  studies 

of  projectile  deceleration  during  hydraulic  ram  have  been  made  by 

2  3 

Stepka  and  Kapple  .   Projectile  position  as  a  function  of  time  was 

measured  using  shadowgraphs.  These  measurements  of  the  velocity  decay 

of  spheres  with  various  diameters,  masses,  and  impact  velocities  indicates 


that  the  energy  deposited  during  shock  formation  is  a  constant  percentage 

of  the  impact  energy  and  is  primarily  a  function  of  the  projectile  impact 

Mach  number  with  respect  to  the  fluid.   Figure  II-2  shows  the  observed 

dependance  of  the  change  in  projectile  velocity  AV  occuring  during  the 

s 

characteristic  shock  formation  time  with  projectile  impact  Mach  number. 
Little  data  was  available  for  times  on  the  order  of  a  few  microseconds, 
hence,  the  data  was  extrapolated  to  small  times  for  these  calculations. 
Impact  Mach  numbers  which  are  less  than  sonic  produce  a  shock  formation 
velocity  decay  which  is  less  than  10  percent  of  the  impact  velocity.  For 
impact  Mach  numbers  greater  than  sonic  the  loss  in  projectile  velocity 
raises  rapidly  to  approximately  80  percent  at  a  Mach  number  of  3.5. 
Figure  II-2  also  shows  the  corresponding  shock  phase  energy  input 
percentage  as  a  function  of  Mach  number.  This  information  was  used  in 
conjunction  with  the  theory  of  Yurkovich  to  estimate  shock  phase  pressures 
during  hydraulic  ram.   These  data  indicate  that  the  shock  phase  becomes 
increasingly  important  as  projectile  impact  Mach  number  increases  since 
a  supersonic  projectile  will  generate  stronger  shock  phase  pressures  than 

a  subsonic  projectile  of  the  same  impact  energy. 

5 
Computer  calculations  based  upon  the  Yurkovich  theory  show  that 

the  shock  is  attenuated  by  spatial  and  viscous  effects  very  quickly  and 

becomes  acoustic.  Figure  II-3  shows  that  the  distance  traveled  into  the 

fluid  and  the  time  from  impact  at  which  the  wave  speed  is  sonic  for  tanks 

filled  with  water  and  with  aviation  fuel.  Figure  II-1+  shows  a  typical 

pressure  distribution  resulting  from  the  impact  of  a  h5   grain  .222  caliber 

small  arms  projectile  with  a  water  filled  tank.  The  wave  produced  becomes 

acoustic  at  approximately  Ik  microseconds  after  impact  and  at  a  distance 

of  one  inch  from  the  impact  point.   The  pressure  distribution  for  times 

larger  than  that  required  for  the  wave  to  reach  sonic  speed  was 

5 


calculated  using  acoustic  wave  theory.   Intense  pressures  are  predicted 
near  the  entry  point  but  the  pulse  is  quickly  reduced  in  magnitude  to 
several  thousand  pounds  per  square  inch  at  a  distance  greater  than  an 
inch.   Figures  II-5  through  7  show  similar  predicted  pressure  distribu- 
tions resulting  from  the  impact  of  a  50  caliber,  12.7  mm  and  lk.5  ram 
threat  in  a  fuel  filled  tank. 
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FIGURE  n-3    RADIUS    AND    TIME  FOR    ACOUSTIC 
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III.  Drag  Phase  Pressures 

Since  the  shock  phase,  in  most  cases  of  interest,  does  not  dissipate 
all  the  available  energy  the  projectile  travels  through  the  fuel  cell 
fluid  until  the  remaining  energy  is  dissipated  or  until  the  exit  wall 
is  impacted.  During  this  drag  phase  of  hydraulic  ram  a  cylindrical 
cavity,  whose  size  is  dependant  upon  the  rate  of  kinetic  energy  loss  along 
the  trajectory,  is  formed.   Projectile  spin  and  angle  of  attack  combine 
to  produce  Magnus  sideforces  which  pertrub  the  projectile  from  a  straight 
path.  Lift  production  usually  causes  projectile  instability  and  the 
projectile  tumbles  a  few  inches  from  the  entry  point.  Large  stagnation 
pressures  in  the  flow  can  cause  projectile  deformation  and  in  some  cases 
the  projectile  does  not  remain  intact.  These  factors  made  an  analytical 
estimation  of  the  drag  phases  pressures  difficult.   Lundstrom  has  analyzed 
the  drag  phase  by  assuming  a  simple  drag  equation  which  results  in  an 
energy  loss  rate  given  by: 

—  =  mBV  /xTX  n  x 

dx  (III-l) 

where  B  is  the  velocity  decay  coefficient  defined  as: 

PfCdA 
2m  (HI-2) 

The  wave  equation  was  used  to  calculate  the  flow  field  resulting  from 
projectile  and  cavity  motion.   The  projectile  and  cavity  were  approximated 
by  a  distribution  of  sources  along  the  trajectory.  The  velocity  potential 
was  expressed  as  an  integral  function  of  the  source  strength  distribution, 
spatial  position,  and  time.   The  source  strength  distribution  was  estimated 
by  equating  the  drag  phase  energy  deposited  by  the  projectile  to  the  sum 
of  the  work  done  by  the  difference  in  ambient  and  cavity  pressure  and  the 
kinetic  energy  contained  in  the  tank  fluid.   Restraining  effects  of  the 
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fuel  cell  walls  on  cavity  growth  are  neglected  and  this  method  cannot 
predict  source  strength  during  cavity  collapse.   The  pressure  fluid 
resulting  from  the  source  distribution  was  calculated  from  the  change  in 
potential  with  time.  Wave  reflections  from  fuel  cell  walls  were  calculated 
using  the  method  of  images  but  waves  reflected  from  the  cavity  surface  are 
ignored.   Fuel  cell  pressure  predictions  are  very  dependant  upon  the 
assumed  motion  of  the  tank  walls.  For  times  less  than  several  hundred 
microseconds  (which  is  characteristic  of  the  duration  of  the  shock  and 
drag  phase)  entry  wall  reflected  waves  are  of  primary  importance. 
Lundstrom's  analysis  can  predict  fuel  cell  pressures  for  two  types  of 
wall  response:  the  rigid  or  stationary  wall  and  the  free  surface.  Thin 

walled  fuel  cell  internal  pressures  may  be  estimated  by  assuming  the  walls 

7 

are  free  surfaces.  Cole  has  shown  that  an  exponential  plane  wave 

impinging  upon  thin  walls  that  are  typical  of  aircraft  fuel  cell  construc- 
tion are  reflected  very  nearly  as  if  the  wall  were  a  free  surface.   This 
assumption,  however,  eliminates  accurate  prediction  of  wall  pressure 
loadings  that  are  necessary  to  predict  wall  motion  and  stress  during 
hydraulic  ram. 

Accurate  prediction  of  fluid  pressures  near  the  walls  requires  the 
solution  of  a  coupled  fluid- structure  interaction  problem.  Ball's 
application  of  piston  theory  to  this  problem,  however,  allows  the  structure 
and  fluid  equations  to  be  uncoupled.   This  theory  requires  only  the 
knowledge  of  the  pressure  and  velocity  that  would  exist  if  the  wall  were 
not  there,  including  earlier  reflections  from  other  walls  and  free 
surfaces.  This  incident  pressure  and  velocity  information  can  be  obtained 
from  Lundstrom's  analysis  and  wall  strains  during  hydraulic  ram  may  be 
estimated. 
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Experimental  Results 

3  11  12 
A  series  of  experiments  '   '  "  were  undertaken  using  a  simple  tank 

configuration  to  determine  hydraulic  ram  pressure  loadings  for  two 

.222  caliber  threat  energy  levels.  The  test  fuel  cell  was  cubical  with 

inside  dimensions  equal  to  seventeen  inches,  and  was  constructed  by 

welding  one  quarter  inch  by  three  inch  aluminum  angle  sections  to  form  a 

frame.  Side  and  bottom  walls  were  one  inch  thick  plexiglass  and  the  entry 

and  exit  walls  were  bolted  to  the  frame  so  that  various  thicknesses  and 

materials  could  be  tested.  The  top  was  open  and  bullet  entry  into  the 

tank  with  minimum  loss  of  energy  was  accomplished  by  providing  a  one 

inch  dimater  hole  in  the  center  of  the  entry  wall.  A  .222  caliber  rifle 

was  used  to  fire  copper- jacketed  lead  projectiles  whose  mass  and  impact 

velocity  was  adjusted  to  obtain  the  7,^+93  in- lb  and  12,323  in- lb  impact 

energy  levels  tested.   Spark  shadowgraphs  were  taken  during  hydraulic  ram 

to  determine  shock  phase  wave  shape,  shock  position  and  projectile  location 

as  a  function  of  time.  Shock  phase  wave  shape  was  found  to  be  hemispherical 

and  the  shock  radius  was  proportional  to  the  0.9  power  of  time  for  the 

projectiles  studied.  Projectile  position  data  was  used  to  evaluate  the 

projectile  velocity  decay  coefficient  B  used  in  the  Lundstrom  drag  phase 

analysis.  These  data  were  used  to  predict  shock  and  drag  phase  internal 

tank  pressures.  Kistler  603H  quartz  pressure  transducers  were  mounted  in 

the  tank  to  measure  internal  pressure  at  a  given  spatial  position  verses 

time.  Pressures  were  measured  along  a  line  inclined  50  from  the  vertical 

at  radii  varying  from  two  to  eight  inches  from  the  impact  point. 

Holm   measured  the  internal  pressures  resulting  when  a  0.5  inch 

thick  steel  plate  was  used  as  the  entry  wall.  This  stimulates  the  case 

of  a  rigid  entry  wall.  Figures  III-l  through  3  show  a  comparison  of  the 
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estimated  and  measured  pressure  for  the  low  energy  level  threat  and  Figure 
111-4  through  6  show  the  higher  energy  level  threat.   The  calculations 
predict  two  pressure  pulses  of  very  different  time  scales.  The  shock 
phase  pulse  is  large  in  magnitude  but  very  short  in  duration  (less  than 
25  microseconds).  The  drag  phase  pulse,  however,  has  approximately 
one-half  peak  shock  phase  magnitude  and  a  duration  of  approximately 
one  hundred  microseconds.  The  experimental  data  does  not  exhibit  the 
characteristics  of  the  predicted  shock  phase  pressures.   The  transducer 
used  has  a  rise  time  of  one  microsecond  but  the  charge  amplifier  used  to 
condition  the  transducer  output  has  a  frequency  response  of  300  KH  .  It 
was  expected,  therefore,  that  the  shock  phase  pulse  would  not  be  accurately 
reproduced  at  the  output.  Yeteerbium  piezoresistive  transducers  with 

response  characteristics  suitable  for  measurement  of  shock  phase  pulses 

13 
have  been  developed  at  the  Stanford  Research  Institute   but  these  thin 

film  gages  are  destroyed  in  a  few  microseconds  after  shock  front  passage. 

Until  better  pressure  transducers  systems  with  rise  times  of  approximately 

10  nanoseconds  are  developed  and  become  commercially  available  the  validity 

of  shock  phase  pressure  predictions  cannot  be  established.  The  calculation 

of  shock  phase  pressure  assumes  a  rigid  entry  wall.  Wall  motion  in 

response  to  the  shock  phase  loading  may  produce  rarefraction  waves  which 

significantly  distort  the  predicted  distribution  of  pressure  near  the  wall. 

The  experiments  do  exhibit  characteristics  similar  to  the  predicted 

drag  phase  pulse.  The  time  scale  and  peak  pressure  magnitude  of  this 

pulse  is  accurately  predicted  but  the  pulse  shape  is  not.  The  drag  phase 

peak  magnitude  comparison  is  shown  in  Figure  III-7.   The  pulse  shape 

discrepancy  was  due  to  the  assumption  of  a  constant  drag  coefficient 

along  the  trajectory.   The  shadowgraph  records  of  projectile  position 

could  not  be  used  to  determine  projectile  tumbling  behavior  because  the 

17 
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cavity  shadow  obliterated  projectile  attitude  information.   Lundstrom 
has  shown  that  if  detailed  projectile  tumbling  behavior  is  included  the 
pulse  shape  during  the  drag  phase  can  be  matched  with  experimental  results, 
The  present  results  show  that  peak  drag  phase  pressures  may  be  estimated 
accurately  using  only  an  estimate  of  the  average  drag  coefficient  of  the 

projectile  as  it  traverses  the  tank. 

12 
Patterson   has  measured  the  internal  pressures  for  the  same  tank 

but  with  an  .05  inch  thick  aluminum  entry  wall.  Figures  II- 8  and  9 

show  a  comparison  of  the  Lundstrom  prediction  (assuming  the  entry  wall  is 

a  free  surface)  and  these  experiments.  Figure  111-10  shows  that  the  peak 

pressure  versus  radius  is  again  accurately  predicted  by  the  Lundstrom 

analysis. 
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IV.   Fuel  Cell  Entry  Wall  Response  During  the  Shock  and  Drag  Phases 

The  estimation  of  fluid- filled  tank  "wall  response  during  hydraulic 
ram  is  extremely  complex  because  the  structural  and  fluid  equations  of 
motion  are  coupled  by  the  presence  of  a  common  boundary.  An  exact 
prediction  of  wall  strain  and  motion  therefore,  required  the  simultaneous 
solution  of  the  structural  and  fluid  differential  equations  of  motion. 

To  make  the  prediction  of  wall  response  more  tractable,  Ball's  '  ' 

application  of  piston  theory  was  used  to  provide  an  estimate  of  the  fluid 

and  structural  interaction  during  hydraulic  ram.   Piston  theory  provides 

the  correct  solution  to  the  one-dimensional  propagation  of  stresses  in 

an  acoustic  medium  due  to  a  moving  boundary.  This  theory  predicts  that 

the  wall  pressure  loading  p  may  be  calculated  by: 

w 

Pw  =  Pi  +  Pc  (v.-w)  (IV_1} 

where  p.  and  v.  are  the  incident  pressure  and  velocity  of  the  fluid  if 
the  wall  were  not  there,  p  is  the  fluid  density,  c  is  the  speed  of  sound 
in  the  fluid  and  w  is  the  wall  velocity.   The  first  two  terms  in  the  wall 
pressure  loading  are  independant  of  wall  motion,  hence  the  application  of 
piston  theory  uncouples  the  structural  and  fluid  equations  of  motion. 
Wall  response  can  be  computed  using  only  the  structural  equations  modified 
to  include  the  third  term  of  equation  IV-1  if  a  known  prediction  of 
incident  pressure  and  velocity  is  used  as  the  loading.   The  fluid  dynamic 
analyses  of  the  previous  two  sections  was  used  to  provide  estimates  of 
the  incident  fluid  pressure  and  velocity  generated  during  the  shock  and 
drag  phases  of  hydraulic  ram.   The  computer  program  SATANS  was  modified 
to  include  piston  theory  and  was  used  to  calculate  the  resulting  wall 
response  due  to  the  time  varying  loading. 
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Experimental  Results 
A  series  of  experiments   were  undertaken  to  determine  the  entry 
wall  response  to  shock  and  drag  phase  loadings  for  two  threat  energy- 
levels.  7075-T6  Aluminum  plates  of  thickness  0.05,  0.09  and  0.l6  inches 
were  tested.   Each  entry  wall  had  four  EA-13  type  175  ohm  strain  gages 
mounted  2  inches  from  the  impact  point.   One  pair  of  gages  was  oriented 
radially  with  one  gage  on  the  inside  and  the  other  on  the  outside  of  the 
plate.  A  second  pair  of  gages  was  oriented  circumferentially.   The  simple 
constant  current  strain  gage  circuit  shown  in  Figure  IV-1  was  used  to 
measure  the  dynamic  response.   The  recording  system  consisted  of  a 
dual-beam  oscilloscope  with  two  1A7  pre-amp  plug-in  units  with  differential 
input  capability  and  a  Polaroid  scope  camera.  The  effects  of  hydrostatic 
pressure  on  the  wet  gage  output  was  neglected  because  it  was  found  to  be 
only  2  to  5  micro-inches  per  inch  of  indicated  compressive  strain  per 
1000  psi.  The  constant  current  circuit  used  in  this  experiment  has  a 
calibration  equation  given  by 

e"   g  1-(1-T1)«K  (IV-1) 

where  the  efficiency  T)  is  given  by: 

„    R2 

11  = 


VRg  (IV-2) 

where  K  is  the  gage  factor,  I  the  current,  R  the  load  resistance,  and 
R  the  active  gage  resistance.   For  the  circuit  elements  shown  in  Figure 
IV-1  the  calibration  equation  was  simplified  to  the  linear  relationship: 

e  =  0.3529e  (IV-3) 

which  was  used  to  calculate  surface  strain  from  the  measured  voltage  output, 

Two  .222  caliber  projectile  energy  levels  were  tested  for  each  of 
the  three  test  plates.   Shots  taken  at  the  higher  energy  level  resulted 


30 


in  considerable  tank  leakage  and  strain  gage  failure.   High  energy  shots 
with  the  0.05  inch  thick  plate  resulted  in  permanent  plate  deformation 
that  was  not  evident  at  the  lower  energy  level.   Significant  deterioration 
of  the  strain  gage  bond  accompanied  each  shot  at  the  high  energy  level. 
As  a  result,  less  data  are  available  at  that  energy  level. 

Figure  IV- 2  shows  typical  experimental  entry  wall  strain 
measurements  versus  time  for  a  0.09  inch  thick  wall  when  impacted  by  the 
71+93  in- lb  energy  projectile.  Additional  data  and  analysis  are  contained 
in  reference  l6.  The  upper  photograph  shows  that  the  radial  strain  is 
due  to  pure  bending  for  approximately  60  p,  seconds  after  impact.  After 
this  time  some  radial  stretching  is  evident  but  both  traces  return  to 
neutral  at  approximately  200  jj,  seconds.  The  dry  (outside  gage)  continues 

below  the  neutral  reading  because  of  stretching  but  the  relatively  symmet- 
rical variations  from  that  point  on  show  the  continuation  of  back  and 

forth  bending.  The  lower  photograph  shows  that  circumferential  strains 
are  also  primarily  due  to  pure  bending  of  the  plate  for  the  first  60  u 
seconds.   The  circumferential  strains,  however,  do  not  return  to  neutral 
which  is  indicative  of  stretching.   The  maximum  strain  amplitude  of  both 
radial  and  circumferential  traces  was  approximately  l.k   X  10   in/in. 

In  order  to  evaluate  the  adequacy  of  piston  theory  for  predicting 
the  entry  wall  structural  response,  three  different  cases  were  calculated 
using  SATMS.  Each  case  used  an  incident  shock  phase  pressure  loading 
that  was  calculated  using  the  computer  program  discussed  in  section  II 
and  a  drag  phase  pressures  loading  predicted  by  the  computer  program  of 
Lundstrom  discussed  in  section  III. 

The  first  case  assumed  no  fluid-wall  interaction.   The  wall  loading 
was  taken  to  be  the  incident  pressure  p.  and  pc(v.-w)  terms  of  piston 
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theory  were  neglected.   Figures  IV- 3  and  k   show  a  comparison  of  this 
result  with  the  experiment.  Predicted  wall  strains  are  very  large 
compared  with  those  measure  experimentally  because  the  reduction  of  fluid 
pressure  at  the  wall  due  to  wall  motion  has  been  ignored.  The  strain 
predictions  indicate  that  the  entry  wall  should  sustain  a  permanent  set. 
Experimental  entry  walls  show  no  such  behavior.   It  is  clear  that  the 
inclusion  of  the  fluid-wall  interaction  is  necessary  of  reasonable  wall 
response  predictions  are  to  be  made  during  hydraulic  ram. 

The  second  entry  wall  response  assumed  that  the  fluid- wall  interaction 
could  be  modeled  by  piston  theory.  The  wall  loading  was  calculated  by 
p.  +  pcv.  where  the  velocity  v.  was  taken  as  the  normal  component  of  the 
incident  fluid  velocity  at  the  wall.  SATANS  was  modified  by  Ball  to 
include  the  p cw  term  in  the  structural  equation  of  motion.  The  results 
of  these  calculations  are  shown  in  Figures  IV- 5  and  6.  These  figures 
show  that  the  predicted  strains  are  considerably  larger  than  measured. 
Figures  7  through  9  show  the  assumed  loadings  of  the  plate  (p.  and 
p.  +  pcv. )  at  three  times  during  the  response  calculation.  The  net  pressure 
at  the  entry  wall  (p.  +  pc(v.-w))  is  also  plotted  in  these  figures.  An 
examination  of  these  figures  shows  a  very  large  loading  of  the  plate 
center  due  to  the  pcv.  term.   This  large  loading  of  the  plate  center 
is  responsible  for  the  overestimation  of  wall  strains.  The  one- dimensional 
pcv.  term  is  clearly  not  a  good  estimate  of  reflected  wave  momentum  near 
the  plate  center  where  the  flow  is  very  three-dimensional. 

The  third  wall  response  prediction  assumed  that  the  wall  loading 
could  be  given  by  2p. .   This  loading  corresponds  to  assuming  that  the 
pcv.  term  can  be  replaced  by  p.  as  in  the  case  of  waves  reflecting  off  a 
rigid  wall.   Figures  IV-  10  and  11  compare  this  modified  piston  theory 
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prediction  and  the  experimental  results.   Figures  IV- 7  through  9  show 
that  the  large  loading  at  the  center  of  the  plate  has  been  eliminated. 
As  a  result  the  predicted  strains  are  much  closer  to  those  measured. 
Figures  7  through  9  also  show  that  the  modified  piston  theory  net  wall 
pressure  (2p.-pcw)  is  negative  toward  the  center  of  the  plate.   This 
indicates  that  the  wall  has  moved  too  quickly  out  of  the  way  in  response 
to  the  fluid  pressure  loading.   The  resulting  negative  wall  pressure  is 
responsible  for  the  descrepancy  between  measured  and  predicted  strains. 
The  negative  net  pressure  may  be  eliminated  by  limiting  the  wall  velocity 
to  values  for  which  the  net  pressure  is  zero. 

Figures  IV- 10  and  11  show  a  peak  strain  at  35  p,  seconds  which 
corresponds  to  the  arrival  of  the  shock  front  at  the  gage  location.   The 
predicted  strains  decrease  only  to  increase  again  as  the  effects  of  the 
drag  phase  pulse  are  felt.   The  experimental  results  do  not  exhibit  the 
characteristics  of  the  first  calculated  peak  which  indicates  that  the  wall 
does  not  respond  to  the  shock  phase  as  predicted.   The  prediction  under- 
estimates the  wall  strain  due  to  the  drag  phase  because  of  the  existance 
of  negative  net  pressure  at  the  wall.   The  agreement  of  theory  and  experi- 
ment is,  however,  much  better  than  any  previous  theory  has  obtained.   It  is 
concluded,  therefore,  that  the  modified  piston  theory  can  be  used  to  give 
reasonable  estimates  of  fuel  cell  wall  response  during  hydraulic  ram. 
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